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Abstract 



In this paper, we provide algorithms to rank and unrank certain degree-restricted classes 
of Cayley trees. Specifically, we consider classes of trees that have a given degree sequence 
or a given multiset of degrees. Using special properties of a bijection due to Egecioglu and 
Remmel [3] , we show that one can reduce the problem of ranking and unranking these classes 
of degree-restricted trees to corresponding problems of ranking and unranking certain classes 
of set partitions. If the underlying set of trees have n vertices, then the largest ranks involved 
in each case are of order n! so that it takes 0{nlog{n) bits just to write down the ranks. 
Our ranking and unranking algorithms for these degree-restricted classes are as efficient as 
can be expected since we show that they require 0{n^log(n) bit operations if the underlying 
trees have n vertices. 

1 Introduction 

In computational combinatorics, it is important to be able to efficiently rank, unrank, and 
randomly generate (uniformly) basic classes of combinatorial objects. A ranking algorithm for a 
finite set 5 is a bijection from S to the set {O,-- - jIS*! — 1}. An unranking algorithm is the inverse 
of a ranking algorithm. Ranking and unranking techniques are useful for storage and retrieval of 
elements of S. Uniform random generation plays a role in Monte Carlo methods and in search 
algorithms such as hill climbing or genetic algorithms over classes of combinatorial objects. 
Uniform random generation of objects is always possible if one has an unranking algorithm 
since one can generate, uniformly, an integer in {0, • • • ,\S\ — 1} and unrank. 

Given the set V = {!,••• ,n}, we consider the set Cn of trees with vertex set V. These 
trees are sometimes called Cayley trees and can be viewed as the set of spanning trees of the 
complete graph Kn. Ranking and unranking algorithms for the set Cn have been described by 
many authors. Indeed, efficient ranking and unranking algorithms have been given for classes 
of trees and forests that considerably generalize the Cayley trees (e.g., [3], [4], [5], [6], [7]). 

In this paper we consider more refined problem, namely, ranking and unranking subsets of 
Cn with specified degree sequences or a specified multisets of degrees. The resolution of this 
refined problem hinges on having a bijection with certain very special properties. Let Cn,i be 
the set of directed trees on V that are rooted at 1. That is, a directed tree T G Cn,i has all 
its edges directed towards its root 1. We replace Cn with the equivalent set Cn,i- For any tree 
T £ Cn, Y17=i degrii) = 2n — 2. Let s = {si, . . . , Sn) be any sequence of positive integers such 



that X]r=i Si = 2n — 2. Then we let Cn^s = G Cn,i '■ {degxiX), • • • , degxin)) = s}. It will 
easily follow from our results in section 2 that 



^Sl - 1, . . . , s„ - 1 

Similarly if 5 = {1'^^ , • • • , (n — 1)°"-^} is a multiset such that - i = 2n — 2 and Y^a=i — 

n, then we let Cn,s = G (7„,i : {degT{l), • • • , degT{n)} = S}. It is easy to see from (llj) that 

ic„.i=( " )( X 

\ai, . .., an) V^i - 1, • • • , Sn - 1/ 
In [5], Egecioglu and Remmel constructed a fundamental bijection between the class of 
functions Tn = {/ : {2, ... ,n — 1} — )• {!,... ,?^}} and the set C^^i. We shall exploit a key 
property of this bijection which is that for any vertex 1 + equals that degree of i in 

the tree T = 0(/) when is considered as an undirected graph. This property allows us to 
reduce the problem or ranking and unranking trees in C„^^ or Cn^s to the problem of ranking 
and unranking certain set partitions. We then use known techniques for ranking and unranking 
set partitions [Tl]. 

Note that both the |C„_^| and jCn^s'l can be as large as as the order of n! so that the numbers 
involved in ranking and unranking will require on the order of nlog(n) bits just to write down. 
We shall show that our algorithms for the ranking and unranking algorithms for C„^^ and Cn^s 
are as efficient as can be expected in that given a tree T, it will require at most O(n^) comparisons 
of numbers y < n plus 0(n) operations of multiplication, division, addition, substraction and 
comparision on numbers x < |C„^^| (x < ICn^sl) to find the rank of T in Cn^g {Cn,s) so that it 
will take 0{n^log{n)) bit operations for ranking in either Cn,s or Cn,s- Similarly, we will show 
that that the unranking algorithms Cn,s or Cn,s will reqire at most 0{'n?log{n)) bit operations. 

The outline of this paper is as follows. In Section 2, we describe the bijection Q : Fn ^ Cn,i 
of [3] and discuss some of its key properties. In Section 3, we show that both © and can 
be computed in linear time. This result allows us to reduce the problem of efficiently ranking 
and unranking trees in Cn^g or Cn,s to the problem of efficiently ranking and unranking certain 
classes of set partitions. In Section 4, we shall give ranking and unranking algorithms for the 
classes of set partitions corresponding the sets Cn^g and Cn,s- 



2 The G Bijection and its Properties 

In this section, we shall review the bijection © : Cn,i due to Egecioglu and Remmel [3] 

and give some of its properties. 

Let [n] = {1,2,..., n}. For each function / : {2, . . . , n — 1} — [n], we associate a directed 
graph /, graph{f) = {[n],E) by setting E = {{i,f{i)) : i = 2, . . . , n — 1}. Following [13], given 
any directed edge where 1 < i,j < n, we define the weight of (i, j), W{{i,j)), by 

«'«'•■'■»={ S\t'>;' (^) 

where pi, qi, Si,ti are variables for i = 1, . . . , n. We shall call a directed edge (z, j) a descent edge 
'd i > j and an ascent edge if i < j. We then define the weight of any digraph G = {[n],E) by 

W{G)= n W{{i,j)). (4) 

(i,i)6E 
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A moment's thought wiU convince one that, in general, the digraph corresponding to a function 
/ G Tn will consists of 2 root-directed trees rooted at vertices 1 and n respectively, with all edges 
directed toward their roots, plus a number of directed cycles of length > 1. For each vertex v 
on a given cycle, there is possibly a root-directed tree attached to v with v as the root and all 
edges directed toward v. Note the fact that there are trees rooted at vertices 1 and n is due to 
the fact that these elements are not in the domain of /. Thus there can be no directed edges out 
of any of these vertices. We let the weight of /, be the weight of the digraph graph{f) 

associated with /. 

To define the bijection B, we first imagine that the directed graph corresponding to f e T 
is drawn so that 

(a) the trees rooted at n and 1 are drawn on the extreme left and the extreme right respectively 

with their edges directed upwards, 

(b) the cycles are drawn so that their vertices form a directed path on the line between n and 

j, with one back edge above the line, and the root-directed tree attached to any vertex on 
a cycle is drawn below the line between n and 1 with its edges directed upwards, 

(c) each cycle Cj is arranged so that its maximum element rrii is on the right, and 

(d) the cycles are arranged from left to right by decreasing maximal elements. 
Figure 1 pictures a function / drawn according to the rules (a)-(d) where n = 23. 




Figure 1: The digraph of a function 

This given, suppose that the digraph of / is drawn as described above and the cycles of 
/ are ci(/), . . . ,Ca(/), reading from left to right. We let ^^.(j) and lci{f) denote the right and 
left endpoints of the cycle Ci{f) for i = l,...,a. Note that if Ci{f) is a 1-cycle, then we 
let = ^Cj(/) be the element in the 1-cycle. ©(/) is obtained from / by simply deleting 

the back edges {r^.^(j-jj,^.^(^f^) for i = l,...,a and adding the directed edges (rc.(/), ^ci+i(/)) fo^" 
z = 1, . . . , a — 1 plus the directed edges (n, lci{f)) cind (r^^jj), 1). That is, we remove all the back 
edges that are above the line, and then we connect n to the lefthand endpoint of the first cycle, 
the righthand endpoint of each cycle to the lefthand endpoint of the cycle following it, and we 
connect the righthand endpoint of the last cycle to 1. For example, ©(/) is pictured in Figure 2 
for the / given in Figure 1. If there are no cycles in /, then @{f) is simply the result of adding 
the directed edge (n, 1) to the digraph of /. 
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Figure 2: e(/) 



To see that © is a bijection, we shall describe how to define 0~ . The key observation is 
that we need only recover that the directed edges (rc.(j), /ci+i(/)) for i = 1, ■ ■ ■ ,a — 1. However 
it is easy to see that ^(.^(j) = mi is the largest element on the path from n to 1 in the tree &{/)■ 
That is, mi is then largest element in its cycle and by definition, it is larger than all the largest 
elements in any other cycle so that nii must be the largest interior element on the path from n 
to 1. Then by the same reasoning, rc2(/) = m2 is the largest element on the path from mi to 1, 
etc. Thus we can find mi, ... , mt- More formally, given a tree T G Cn,i, consider the path 

mo = n,xi, . . . ,mi,X2, . • .m2, . . . ,xt, . . . ,mt,l 

where m^ is the maximum interior vertex on the path from m^-i to 1, 1 < z < i. If (mj_i,mi) is 
an edge on this path, then it is understood that Xi, . . . ,mi = rrii consists of just one vertex and 
we define Xi = rrii. Note that by definition mo = n > mi > ... > nit- We obtain the digraph 
9~^(r) from T via the following procedure. 

Procedure for computing 0~^(T) : 

(1) First we declare that any edge e of T which is not an edge of the path from n to j is an edge 

of G-i(r). 

(2) Next we remove all edges of the form (mt, 1) or (mi_i,Xi) for 1 < i < t. 
Finally for each i with 1 < i < we consider the subpath Xi, . . . ,mi. 

(3) If rrii = Xi, create a directed loop {rrii, rrii). 

(4) If rrii 7^ Xi, convert the subpath Xi, . . . ,mi into the 
directed cycle Xi, . . . ,mi,Xi. 



Next we consider two important properties of the bijection O. First Q has an important 
weight preserving property. We claim that if Q{f ) = T, then 

qnhWif) = W{T). (5) 

That is, by our conventions, any backedge (r^^(^f),lf..(^f^) are descent edges so that its weight is 
Qr (t\U 1,^- Thus the total weight of the backedges is 

a 
i=l 
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Our argument above shows that aU the new edges that we add are also descent edges so that 
the weight of the new edges is 



a-l 



i=l i=l 

Since ah the remaining edges have the same weight in both the digraph of / and in the digraph 
e(/), it follows that qnhWif) = T^(G(/)) as claimed. 
It is easy to see that 

n-l 

Yl W{f) = Y[[qi{ti + --- + U)+pi{si+i + --- + Sn)]. (8) 

Thus we have the following result which is implicit in [3] and it explicit in |13j . 
Theorem 1. 

n-l 

W{T) = qnti'[[[qi{ti + ---+ti)+pi{si+i + --- + Sn)]. (9) 



rec„ 



i=2 



Next we turn to a second key property of the Q bijection. It is easy to see from Figures 
1 and 2 that deleting the back edges (?'ci{/)) ^Ci(/)) for z = 1, . . . , a in graph{f) and adding the 
directed edges ('^Cj(/)) ^Ci+i(/)) for i = 1, . . . , a — 1 plus the directed edges (n, ^ci(/)) and (rc^(j), 1) 
to get Q{f) does not change the indegree of any vertex except vertex 1. That is, 

indegg^aph{f){i) = indegQ(^f){i) for z = 2, . . . , n. (10) 

It is also easy to see that in going from graph{f) to @{f), the indegree of vertex 1 increases by 
1, i.e., 

1 + indeggraph(f){i) = indegQ^f){l). (11) 

When we consider 0(/) as an undirected graph T, then it is easy to see that degrii) = 
outdegQ(^fj + indegQ(^jy Thus since the outdegree of i in B(/) is 1 if i 7^ 1 and the outde- 
gree of 1 in 0(/) is zero, equations (10) and (11) imply the following Theorem. 

Theorem 2. Suppose that T is the undirected tree corresponding to G(/) where f G J-n, then 
for i = 1, . . . ,n, 

degT{i) = l + \.r'{i)\. (12) 

ProofBy our definition of graph{f), it follows that indegg^aph{f ){''') = for i = 1, . . . , n. 

Thus by ( 10 ) , for z = 2, . . . , n, 

degrii) = outdegQ(^f){i) + indege(^f){i) 
= 1 + mdec/0(/)(i) 
= 1 + indeggraph(f)ii) 

= 1 + ir'wi- 
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Similarly by (11), 



degril) = outdegQ(^f){l) + indegQi^f){l) 
= + indegQ(^f){l) 
= I + indeggraph{f){'^) 

= i + ir'(i)i- 

□ 



3 Construction of the spanning forests from the the function 
table in time 0{n) 

In this section, we shall briefly outline the proof that one can compute the bijections and its 
inverse in linear time. Suppose we are given / G J^n- Our basic data structure for the function 
/ is a list of pairs {i, f{i)) for i = 2, . . . ,n — 1. Our goal is to construct the directed graph of 
Q{f) from our data structure for /, that is, for i = 1, . . . ,n, we want to find the set of pairs, 
{i,ti), such that there is directed edge from i to tj in &{/)■ We shall prove the following. 

Theorem 3. We can compute the bijection Q : J^n Cn,i and its inverse in linear time. 

Proof. We shall not try to give the most efficient algorithm to construct ©(/) from /. In- 
stead, we shall give an outline the basic procedure which shows that one can construct @{f) 
from / in linear time. For ease of presentation, we shall organize our procedure so that it makes 
four linear time passes through the basic data structure for / to produce the data structure for 

e(/). 

Pass 1. Goal: Find, in linear time in n, a set of representatives ti,...,tr of the cycles of 
the directed graph of the function f. 

To help us find ti, . . . ,tr, we shall maintain an array ^[2], yl[3], • • • A[n — 1], where for each i, 
A[i] = {ci,pi, Qi) is a triple of integers such Cj G {0, . . . , n — 1} and {pi, qi} C {—1, 2, • • • , n — 1}. 
The Cj's will help us keep track of what loop we are in relative to the sequence of operations 
described below. Then our idea is to maintain, through the pi and qi, a doubly linked list of the 
locations i in A where Cj = 0, and we obtain pointers to the first and last elements of this doubly 
linked list. It is a standard exercise that these data structures can be maintained in linear time. 

Initially, all the Cj's will be zero. In general, if Cj = 0, then pi will be the largest integer j 
such that 2 < j < i for which Cj = if there is such a j and pi = —1 otherwise. Similarly, we 
set qi > i to be the smallest integer k such that n — 1 > k > i for which = if there is such 
a k and qi = —1 if there is no such k. If C2 > 0, then q2 is the smallest integer j > 2 such that 
Cj = and 52 = —1 if there is no such integer j. If c„_i > 0, then pn-i is the largest integer 
k < ni such that Cfc = and Pn-i = — 1 if there is no such integer k. 

We initialize A by setting A[2] = (0, -1,52), A[i] = (0, i - 1, i + 1) for m + 1 < i < n - 1, 
and A[n — 1] = (0,p„,_i, —1). If 2 < n — 1 then 92 = 3 and Pn-i = n — 2. Otherwise (2 = n — 1), 
these quantities are both —1. 

LOOP(l): Start with ii = 2, setting C2 = 1. Compute /0(2), /i(2), /2(2), . . . , /^i(2), each 
time updating A by setting Cp(2) = 1 and adjusting pointers, until, prior to setting Cjfci(2) = 
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we discover that either 

(1) f^{2) £ {l,n}, in which case we have reached a node in graph{f) which is not in the 
domain of / and we start over again with the 2 replaced by the smahest i for which a = 0, or 

(2) X = f^^{2) already satisfies Cx = 1- This condition indicates that the value x has al- 
ready occurred in the sequence 2, /(2), /^(2), . . . , f^^{2). Then we set ti = f^^{2). 

LOOP(2): Start with 

— Qm+i which is the location of the first i such that Cj — 0, 
and repeat the calculation of LOOPl with 12 instead of ii = 2. In this manner, generate 
/'^(^2), /^(^2), /^(^2)) • • • ) /*'^(^2)5 each time updating A by setting cp^^^-^ = 2 and adjusting 
pointers, until either 

(1) f^^{i2) G {li'i-}, in which case we have reached a node in graph{f) which is not in the 
domain of / and we start over again with the 12 replaced by the smallest i for which a = 0, or 

(2) X = /^^(ia) already satisfies Cx = 2. (This condition indicates that the value x has al- 
ready occurred in the sequence i2, fih), f'^i'i'2), ■ ■ ■ , f^^ (^2)-) Then we set t2 = f^^ [12) ■ 

We continue this process until 52 = —1- At this point, we will have generated ti, . . . , t,-, where the 
last loop was LOOP(r). The array A will be such that, for all2<i<n — l,l<Cj<r identifies 
the LOOP in which that particular domain value i occurred in our computation described above. 

Pass 2. Goal: For i = 1, . . . ,r, find the largest element rm in the cycle determined by ti. 
It is easy to see that this computation can be done in linear time by one pass through the array 
A computed in Pass 1 above. At the end of Pass 2, we set U = f{mi). Thus when we draw the 
cycle containing ti according to our definition of Qj{f), rrii will be right most element in the 
and li will be the left most element of the cycle containing ti. However, at this point, we have 
not ordered the cycles appropriately. This ordering will be done in the next pass. 

Pass 3. Goal: Sort (/i, m-i), . . . , (/fc, ?Tifc) so that they are appropriately ordered according the 
criterion for the bijection 6(/) as described in by condition (a) -(d). 

Since we order the cycles from left to right according to decreasing maximal elements, it is 
then easy to see that our desired ordering can be constructed via a lexicographic bucket sort. 
(See Williamson's book ^14j for details on the fact that a lexicographic bucket sort can be carried 
out in linear time.) 

Pass 4. Goal: Construct the digraph of Q{f) from the digraph of f. 

We modify the table for / to produce the table for ©(/) as follows. Assume that {li, mi), . . . , (Ik^mf^) 
is the sorted list coming out of Pass 3. Then we modify the table for / so that we add entries for 
the directed edges (n, ^i) and (m^, 1) and modify entries of the pairs starting with mi, . . . , ruk 
so that their corresponding second elements are I2, . ■ . ,lk,^ respectively. This can be done in 
linear time using our data structures. 

Next, consider the problem of computing the inverse of 0. Suppose that we are given the 
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data structure of the tree T G Ci, i.e. we are given a set of pairs , {i,ti), such that there is a 
directed edge from i to U in T. Recall that the computation of Q~^{T) consists of two basic steps. 

Step 1. Given a tree T G Cn,i, consider the path 

mo = n,xi, . . . ,mi,X2, . ..m2, . . . ,xt, . . . ,mt,l 

where rrii is the maximum interior vertex on the path from mj_i to 1, 1 < i < t. If (mi_i, nii) 
is an edge on this path, then it is understood that Xi, . . . ,mi = rrii consists of just one vertex 
and we define Xi = rrii. Note that by definition mo = n > mi > . . . > mt- 

First it is easy to see that by making one pass through the data structure for F, we can construct 
the directed path n — > ai — > . . . — > where 1 = o,,.. In fact, wc can construct a doubly linked 
list (n, ai, . . . , Ur-i, 1) with pointers to the first and last elements in linear time. If we traverse 
the list in reverse order, (1, a^-i, . . . , ai, n), then it easy to see that mt = flr-i) "mti is the next 
element in the list (0^-2, ■ ■ ■ ,0-1) which is greater than mt and, in general, having found m^ = ag, 
then mj_i is the first element in the list (a^-i, . . . , ai) which is greater than rrij. Thus it is not 
difficult to see that we can use our doubly linked list to produce the factorization 

mo = n,xi, . . . ,mi,X2, . . .m2, ...,xt,.. .,mt, 1 

in linear time. 

Step 2. We obtain the digraph Q~^{T) from T via the following procedure. 
Procedure for computing QJ^{F) : 

(1) First we declare that any edge e of T which is not an edge of the path from n to 1 is an edge 

of e-i(r). 

(2) Next we remove all edges of the form (m^, 1) or (mj_i, Xj) for 1 < i < i. 
Finally for each i with 1 < i < t, we consider the subpath Xi,...,mi. 

(3) If mi = Xi, create a directed loop {mi, mi). 

(4) If mi Xi, then , convert the subpath Xi,...,mi into the 
directed cycle Xi, . . . ,mi,Xi. 

Again it is easy to see that we can use the data structure for T, our doubly linked list, and 
our path factorization, mo = n,x\, . . . ,m\,X2, . ■ ■m2, ■ ■ ■ ,xt, . ■ ■ ,mt,j to construct the data 
structure for graph{f) where / = 0~^(T) in linear time. □. 

Given that we can carry out the bijection Q and its inverses in linear time, it follows that 
in linear time, we can reduce the problem of constructing ranking and unranking algorithms 
for Cn,i to the problem of constructing ranking and unranking algorithms for the corresponding 
function class J-n- In the next section, we will construct our desired ranking and unranking 
algorithms for the function classes corresponding to sets of trees Cn,g and Cn,s described in the 
introduction. 
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4 Ranking and Unranking Algorithms for Trees with a Fixed 
Degree Sequence. 

Recall that if T S Cn, then Yll=i degrii) = 2n — 2. Let s = {si, . . . , s„) be any sequence of posi- 
tive integers such that Y17=i ~ 2n— 2. Then we define Cn^g = {T £ Cn,i '■ {degT{l), • • • , degrin)) = 
s}. Similarly if 5" = • • ■ , (n — 1)""^^} is a multiset such that Y17=i oti ■ i = 2n — 2 and 
^11=1 Q^i = then we define Cn,s = {T £ Cn,i '■ {degril), • • • , degxin)} = S}. The main goal 
of this section is to construct 'n?log{n) time algorithms for ranking and unranking trees in Cn^g 
and Cn,s- 

So assume that s = {si, . . . , s„) is a sequence of positive integers such that Y17=i Si = 2n — 2. 
By Theorem[2| it follows that if Q{f) = T, then degxii) = 1 + 1/^^(01 i = 1, . . .n. It follows 
that 

Q-HCn,s) = {/ G -^n : (ir'(l)l, . ..,\f-Hn)\) = {si-l,...,Sn- 1)}. (13) 

Since a function / S J-n is clearly determined by the sequence (/~^(1), . . . ,/~-^(n)), it follows 
from our results in Section 2 that the problem of finding an algorithms to rank and unrank 
trees Cn^g can be reduced to the problem of ranking and unranking ordered set partitions 
(vTi, . . . , TTn) of {2, . . . , n — 1} where the sizes of the sets are specified. That is, we need to 
find an algorithm to rank and unrank ordered set partitions in n„^, the set of all sequences 
of pairwise disjoint sets (vri, . . . , 7r„) such that Ufc=i '''"A: — — 1} and Itt^I =5^ — 1 

for k = 1, . . . ,n. The total number of elements in n„^^ is clearly the multinomial coefficient 

L-CL-i) = (r-i)r'.i""'^) •••(""'"& ^'~'^)- Thus our first step is to develop a sim- 
pie algorithm to rank and unrank objects corresponding to a product of binomial coefficients 

n-=iao- 

For a single binomial coefficient (^) , we shall rank and unrank the set VFn^k of decreasing 
functions / : {1,...,A;} — t- {l,...,n} relative to lexicographic order. A number of authors 
have developed ranking and unranking algorithms for T>J^n,k- We shall follow the method of 
Williamson |14j . First, we identify a function / : {1, . . . ,k} — >• {1, . . . , n} with the decreasing 
sequence (/(I), . . • , f{k)) where n > /(I) > . . . > f{k) > 1. We can then think of the sequences 
as specifying a node in a planar tree Tx>j^„ ^ which can be constructed recursively as follows. 
At level 1, the nodes of Tjyjr^^ ^ are labeled k, . . . ,n from left to right specifying the choices for 
/(I). Next below a node j at level one, we attach a tree corresponding to T-dj^ _^ where a 
tree Tx>Ta i consists of a tree with a single vertex labeled a. Figure 3 pictures the tree T^j-g 3 . 




Figure 3: The tree T^jr^ 3 

Then the decreasing sequence (6,2,1) corresponds to the node which is specified with an 
arrow. It is clear that the sequences corresponding to the nodes at the bottom of the tree T-j^jr^ 3 
appear in lexicographic order from left to right. Thus the rank of any sequence {f{l),...,f{k)) S 



9 



VJ-n^k is the number of nodes at the bottom of the tree to the left of the node corresponding to 
(/(I), . . . , f{k)). Hence the sequence (6, 2, 1) has rank 10 in the tree Tx>j^g ^. 

This given, suppose we are given a sequence (/(I), • • • , f{k)) in Tx>j^^ ^. Then the number of 

leaves in the subtrees corresponding the nodes k, . . . , /(I)— 1 are respectively (^ij) , (^^i) , • • • , {^^^1 
Thus the total number of leaves in those subtrees is 

k-lj \k-lj \ k-l J V k 

Here we have used the well known identity that Yl^s=X~i ik-i) ~ (fc)' follows that the rank of 
(/(I), . . . , f{k)) in TvT^ , equals (^(1)-^) plus the rank of (/(2), . . . , f{k)) in T^,^,,,,., . The 
following result, stated in [14J, then easily follows by induction. 

Theorem 4. Let f : {1, . . . ,k} — t- {1, . . . ,n} be a descreasing function. Then the rank of f 
relative to the lexicographic order on 'DJ^n,k is 

It is then easy to see from Theorem |4j that the following procedure, as described by 
Williamson in [TJ], gives the unranking procedure for DJ^n^k- 

Theorem 5. The following procedure UNRANK{m) computes 

f = if {!),..., f{k)) 
such that Rank-jyT^ ^{f) = m for any I < k < n and < m < (I) - I. 

Procedure UNRANK(m) 

initiaUze m' := m, t := 1, s := k; (1 < k < n, < m < (^') — Ij 
while t < k do 



begin /(t) — 1 = max{y : (^) < m'}; 
m' :=m'- p)-l); 
^ := t + 1; 
s := s — 1; 

end 

Note that {VTa^^^bi ^ '^-^02,^2 x ••• x T^J^at,bt\ ~ Y\l=i iV)- Thus we can use VJ^aiM ^ 
T>J^a2,b2 X • • • X T^^atM °^ objccts Corresponding to a product of binomial coefficients. 

We shall idenitfy an element 

(/i, • • • , /<) G T^^a^M X ^-^02,62 X • • • X VFatM 

with a sequence 

(/i(l), . . . , /2(1), . . . , 72(62), . . . , /i(l), . . . , ftih)) 
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and rank these sequences according to lexicographic order. To define our ranking and unranking 
proceedure for this set of sequences, we first need to define a product relation on planar trees. 

Given a rooted planar tree T, let L{T) be the numbers of leaves of T and Path{T) be the 
set of paths which go from the root to a leaf. Then for any path p G Path{T), we define the 
rank of p relative to T, rankT{p), to be the number of leaves of T that lie to the left of p. 

Given two rooted planar trees Ti and T2 , Ti<S> T2 is the tree that results from Ti by replacing 
each leaf of Ti by a copy of T2, see Figure 3. If the vertices of Ti and T2 are labeled, then 
we shall label the vertices of Ti (g) T2 according to the convention that each vertex v in Ti 
have the same label in Ti (g) T2 that it has in Ti and each vertex w in a copy of T2 that is 
decendent from a leaf labeled I in Ti has a label {l,s) where s is the label of w in T2. Given 
rooted planar trees Ti, . . . ,Tk where k > 3, we can define Ti T2 • • • by induction as 
(Ti (g) • • • (g) Tk_i) (g Tjfc. Similarly if Ti . . . , Tj. are labeled rooted planar trees, we can define the 
labeling of Ti (g r2 (g • • • (g by the same inductive process. 



Now suppose that we are given two rooted planar trees Ti and T2 and suppose that pi G 
Path{Ti) and p2 G Path{T2). Then wc define the path pi ® p2 in Ti (g T2 which follows pi to 
its leaf I in Ti and then follows p2 in the copy of T2 that sits below leaf / to a leaf I') in 
Ti (g) T2. Similarly, given paths G Tj for i = 1, . . . A;, we can define a path p = pi ® ■ ■ ■ ® pk ^ 
Path{Ti (g T2 (g ■ • • (g Tfe) by induction as (pi (g • • • (g Pk-i) 'Si Pk- 

Next we give two simple lemmas that tell us how to rank and unrank the set of paths in 
such trees. 

Lemma 6. Suppose that Ti, . . . , are rooted planar trees and T = Ti (g r2 (g ■ ■ • (g T/t- Then 
for any path p = pi (g • • • (g pfc G Path{T), 



Proof. We proceed by induction on k. Let us assume that Ti, . . . are labeled rooted planar 
trees. 

First suppose that k = 2 and that pi is a path that goes from the root of Ti to a leaf labeled 
li and p2 goes from the root of T2 to a leaf labeled l2- Thus pi (g) p2 goes from the root of 
Ti (g) T2 to the leaf in Ti and then proceeds to the leaf (li, I2) in Ti 0X2. Now for each leaf 
I' to the left of h in r2, there are L{T2) leaves of Ti (g T2 that lie to left of {h, I2) coming from 





Figure 4: The operation Ti (g r2. 



k 



k 




(15) 
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the leaves of the copy of T2 that sits below /'. Thus there are a total of L(T2) ■ rankj-iipi) such 
leaves. The only other leaves of Ti T2 that lie to left of pi (8> P2 are the leaves of the form 
where I" is to left of p2 in T2. There are rankT2iP2) such leaves. Thus there are a total 
of rankT2{P2) + -^(^2) • rankT^{p\) leaves to left of pi ® p2 and hence 



rankTi(^T2iPi ® P2) = rankT2{P2) + L{T2) ■ rankT^{pi] 



as desired. 

Next assume that (15) holds for k <n and that n > 3. Then 
ranA;Ti®-cg)T„(pi <^ ■ ■ ■ <^ Pn) = 

™"-Vi®-®T„_i)®T„(Pl ■ ■ -(^Pn-l) ®Pn) = 
n— 1 n— 1 

rankx^iPn) + L{Tn)C^rankT^{pj) L(Ti)) = 
j=i i=j+i 

n n 

^rankr^iPj) L{Ti). 
j=i i=j+i 

□. 

This given, it is easy to develop an algorithm for unranking in a product of trees. The proof 
of this lemma can be found in 



Lemma 7. Suppose that Ti, . . . , T/^ are rooted planar trees and T = Ti ® T2 ® ■ ■ ■ ® T^. Then 
given ap £ Path[T) such that rankxip) = r^, p = pi®- ■ -^pt G Path(T) where rankxiipi) = Qi 
and 



ro 



qi Jl L{Ti) + n where < n < Yli=2 HTi), (16) 
1=2 

k 

ri = (?2 n ^(^') + ^2 where < r2 < Jits ^Ti), (17) 
1=3 

'■■ (18) 

rk-2 = qk-iL{Tk) + r^-i where < r^-i < L{Tk) and (19) 

rfc_i = qk (20) 

It follows that ranking and unranking our sequences 

(/i(l), . . . , /i(6i), /2(1), . . . , 72(^2), . . . , /t(l), . . . , ft{ht)) 

coresponding to an element {fi, . . . , ft) G VJ^aiM ^ ^-^a2,fe2 x • • • x T^J^atM^ need only rank 
and unrank the leaves with respect to the tree 

That is, consider a path p = pi P2 ® ■ ■ ■ ® Pt G r(ai,6i),...,(at,6t)- For each i, pi corresponds to a 
sequence (/i(l), . . . , /j(6j)) in VFai,bi and hence p corresponds to the sequence 

(/l(l),...,/l(6l),/2(l),...,/2(62),...,/t(l),...,/i(6i)). 
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We are now in position to give an algorithm to rank and unrank ordered set partitions in n„^^, 
the set of ah sequences of pairwise disjoint sets (vri, . . . , 7r„) such that |Jfc=i '^k = i"^, ■ ■ ■ ,n — 1} 
and IvTfcl = Sfc — 1 for k = 1, . . . ,n. Since the total number of elements in n„_^ is clearly the 

multinomial coefficient = (,"11) (""^^^-1"^^) • • • (""^"^^f ''"^^) , we shall identify 

an ordered set partition vr = (vri, . . . , VTn) with an element 

ifli ■■ ■ ,fn) e 'DTn-2,si-l X 'DTn~2~{si-l),S2-l X • • • X PJ"„_2-(J]^^^ Si-l),s„-l 

as follows. Suppose that n = 12 and s = (si, . . . , S12) = (1, 1, 3, 1, 4, 1, 3, 1, 2, 1, 3, 1). Note that 
Y^ll^ Si = 2(12) -2 = 22 so that this is a possible degree sequence for a tree in Ci2- For example, 
the tree Tq G Ci2,i pictured in Figure 5 has this degree sequence when considered as a tree. 
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Figure 5: The tree Tq 

Note that in this case, {si - 1, . . . , S12 - 1) = (0,0,2,0,3,0,2,0, 1,0,2,0). Also in Figure 5, 
we have pictured the graph of / = 0~^(To) and in this case 

^To = (r'(l),...,r'(12)) = (0,0,{lO,7},0,{ll,8,4},0,{9,3},0,{6},0,{5,2},0). 

It will be more efficient for our ranking and unranking procedure to order the set partition by 
increasing size of the parts. Thus we will make one pass through the to extract the |/~^(^)| for 
each i and its relative rank for those parts of the same size. In this case, we would produce the 
following list 

((1, 0, 0), (2, 0, 1), (3, 2, 0), (4, 0, 2), (5, 3, 0), (6, 0, 3), (7, 2,1), (8,0,4), (9, 1,0), (10,0, 5), (11, 2, 2), (12, 0, 6)). 

Here for example, then entry (7, 2, 1) means that the size of /~^(7) is 2 and there is one element 
i < 7 such that |/~^(i)| = 2. We can then do a lexicographic bucket sort to produce a list of 
the elements according to lexicographic order on the last two entries of this list in linear time, 
see [13]. Thus in our example we would produce the following list. 

((1, 0, 0), (2, 0, 1), (4, 0, 2), (6, 0, 3), (8, 0, 4), (10, 0, 5), (12, 0, 6), (9, 1, 0), (7, 2, 1), (3, 2, 0), (11, 2, 2), (5, 3, 0)>. 

The set partition corresponding to this order is 



V-i, V^, 



U6}, {10, 7}, {9,3}, {5,2}, {11,8, 4}). 

We can ignore the 0's and just consider the reduced partition 

^To = ({6}, {10, 7}, {9, 3}, {5, 2}, {11, 8, 4}). 

More generally, let D = {2, . . . ,n — 1}. Let vf = {Bi, . . . , Bk) be an ordered set partition of 
D where each block Bi is nonempty and ordered in decreasing order, Bi = {bi^i > . . . > h^t^}, 
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that comes from some element vr in ni2,s as described above. Let (ri^i, . . . , ri^ti) be the ordered 
sequence of ranks of the respective (61,1, . . . in D. In general, let (rj^i, . . . ,rj^t-) be the 

ranks of the respective . . . , 6^,* J in D \ U^~}^Bj. For tttq, 

= (5) can be considered an element of D-Fio,i 
('^2,1, ^^2,2) = (8,5) can be considered an element of 'DTg^2 
('^3,1; '^3,2) = (6,2) can be considered an element of VT-j^2 
(r4^i,r4_2) = (3, 1) can be considered an element of VF^^2 
(^5,1) ^5,2, ^"5,3) = (3, 2, 1) can be considered an element of VF^^z 

Thus we can think of ttto as the sequence (6,8,5,6,2,3,1,3,2,1) coming from an element of 
^-^10,1 X ^-^9,2 X T^^7,2 X T^^h,2 X 25-^3^3 or as a leaf in the tree Tx>r^^^^ (g) Ixij-g ^ ® T-djt^ ^ 
T-dTt, 2 ® '^T'j^s 3 • Note that the size of the trees needed for the product lemma. Lemma [6| are 

l^VTsJ = (3) = 1' 

® TvT,,2 ® T^^^sJ = (2) " (2) " (3) = 210, 
1^25.^9,2 TvT7,2 ® Tv:f5,2 ® Tv:f5,2 ® ^I'^s.al = (2) ■ (2) ■ (2) ■ (3) = "^560, 

iTb^io,! ® TvT,,, TvT,,, TvT,^, TvT,^,\ = (\°) Q ■ Q) ■ H) ■ H) = 75600. 
Thus we can apply Lemma |6] and conclude that the 

rank{TfTQ) 



By Theorem [4j we have that 

ran/coj-io_i((5)) = 
ran/ci5jrg2((8,5)) = 

ran/cx^jr, 2((3, 1)) = 

Thus 

rank{Tfo) = (4 x 7560) + (25 x 210) + (11 x 10) + (1 x 1) = 35, 601. 
Hence the tree Tq pictured in Figure 5 has rank 35,601 among all the trees 

T S (^12,(0,0,2,0,3,0,2,0,1,0,2,0} • 

If we are given, the degree sequence s, we can assume that we preprocess the sizes of the trees 
needed to apply the product lemma. Lemma [6j Thus we need only compute 0{n) products, 
additions, and multinomial coefficients. Again given, s, we can construct a table of all the 
possible binomial coefficients that we need as part of the preprocessing. Thus to find the rank of 
tree To requires only a linear number of muliplications, additions and table look ups for numbers 
X < |Cn,sl- Since each x requires at most 0{nlog{n)) bits, it is easy to see that these operations 
require at most 0{v?log{n)) bit operations. The only other contribution to the complexity of 



= rank-Djr^^ ^{{b)) x 7560 
+ranA;Dj-9 2((8,5)) x 210 
+rarafcDjr7 2((6, 2)) x 10 
+ran/cx)j-5 2((3, 1)) x 1. 

(v)-(^r)-— • 

('2 VCi')='°+^="' 
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the algorithm is the time it takes go from the representation of the tree to the corresponding 
rank sequences 

{riA, ri,ti , r-2,1, . . . , r2,t2,- ■ ■ , rk,i, • • • , rk,tk)- 

It is easy to see from the fact that we can compute in linear time that we can start with a 
tree and produce the ordered set partition ttt = {Bi, . . . , Bk) in linear time. Thus to complete 
our analysis of the the complexity of the ranking prodeedure, we need to know the complexity 
of the transformation between the (Bi, . . . , Bk) and ranks 

• • • , n,ti , r2,i, T2,t2, • • • , rk,i, . . . , rk,tk)- 

Lemma 8. Let D = {2, . . . ,n — 1}. Let {Bi, . . . , B^) he an ordered partition of D where each 
block Bi is nonempty and ordered in decreasing order, Bi = {bi^i, . . . , ^i,*, ). Let {ri^i, . . . , ri^n) be 
the ordered sequence of ranks of the respective . . . in D. In general, let {ri^i, . . . ,ri^ti) 

be the ranks of the respective . . . , in D \ U^jZ}iBj. Given the sequences (rj^i, . . . , Tj^t-), 
i = 1,. . . ,k, the sets partition Bi, . . . ,Bk can be constructed in worst case time 0{n^log{n)). 
Conversly, given the set partition Bi, . . . ,Bk the sequences (rj^i, . . . , rj^tj, i = 1, . . . ,k can be 
constructed in worst case time 0{n?log{n)). 

Proof. First we can make one pass through the list and set bij : bij — 1 so {Bi, . . . ,Bk) 
become an ordered partition of {1, . . . , n — 2}. Conversly, we can go from an ordered set partition 
. . . , Bk) of {1, ... , n— 2} to an ordered set partition of {2, ... , n— 1} by setting 6jj := bi^j + 1. 
Thus there is no loss in assuming that {Bi, . . . ,Bk) is an ordered partition of {1, . . . , n — 2}. 

This given, it will then be the case that . . . , = • • • , T~i,ti)- Then it will take 
0{n) comparisions of numbers less than or equal to n to construct a sequence /(I), . . . ,f{n) 
where /(i) = : bi,j < i,j = Then it will take nlog{n) steps to create 

the sequences . . . , bi^tj for i = 2, . . . , /c where bij = bij — f{bij). It then easily follows 
that we have reduced the problem to finding the transformation from the ranks (rj^i, . . . , rj j.), 
i = 2, . . . , A;, to set partitions B2, ■ ■ ■ , Bj^ which we can do by recursion. 

It then easily follows that given the set partition Bi, . . . , Bk the sequences (rj^i, . . . , rj^^J, 
i = 1, . . . , k can be constructed in worst case time 0{n?'log{n)) and that given the (r^^i, . . . , rj^^.), 
i = 1, . . . , A;, we can construct the set partition B\, . . . ,Bk in worst case time 0{n'^log{n)). □ 

It then follows that our ranking procedure for Cn^g requires 0{v?log{n) bit operations. 

The unranking procedure for comes from simply reversing the ranking procedure using 
Theorem [5] and Lemma [7| Again we will exhibit the procedure by finding the tree Ti whose 
rank is 50,005 in C12, (0,0,2,0,3,0,2,0, 1,0,2,0} • The first step is to carry out the series of quotients and 
remainders according to Lemma [7} In our case, this leads to the following calculations. 

50, 005 = (6 X 7560) + 4645 
4645 = (22 X 210) + 25 
25 = (2x10) + 5 
5 = (5xl)+0. 

It then follows that we can construct the sequence corresponding to tttx by concatonating the 
sequences ui , . . . , where 

1. Ml is the decreasing function of rank 6 in DJ^io,i, 
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2. U2 is the decreasing function of rank 22 in DTg^2, 

3. M3 is the decreasing function of rank 2 in T)Tt^2i 

4. n4 is the decreasing function of rank 5 in T>T^^2 and 

5. M5 is the decreasing function of rank in T)T^^z- 

It is clear that the sequence of rank 6 in VFiq^i is (7). 

To find the element (/(l),/(2)) of rank 22 in T)J-^^2^ we use the procedure in Theorem [sj 
We start by setting m' := 22 and s = 2. Since (^) = 21 < 22 < (^) = 28, then /(I) -1 = 7 
and hence /(I) = 8. Then we set m! := 22 — 21 = 1 and s = 1. Since 1 — (^) = 0, we get that 
/(2) - 1 = 1 or /(2) = 2. Thus (8, 2) has rank 22 in PJ"9,2. 

To find the element (/(I), /(2)) of rank 2 in 'DTt^2^ we again use the procedure in Theorem 
[sj We start by setting m! := 2 and s = 2. Since (2) = 1 < 2 < (2) = 3, then /(I) -1 = 2 
and hence /(I) = 3. Then we set m' := 2 — 1 = 1 and s = 1. Since 1 — (|^) = 0, we get that 
/(2) - 1 = 1 or /(2) = 2. Thus (3, 2) has rank 2 in VT'j^2- 

To find the element (/(I), /(2)) of rank 5 in T)T^^2^ we again use the procedure in Theorem 
[sj We start by setting m! := 2 and s = 2. Since (2) = 3 < 5 < (2) = 6, /(I) — 1 = 3 and hence 
/(I) = 4. Then we set m' := 5 - 3 = 2 and s = 1. Since 2 - (^) = 0, we get that /(2) -1 = 2 
or /(2) = 3. Thus (4, 3) has rank 2 in VT^^2- 

Finally there is only one element in T)T^ 'i which is (3,2,1). Since the last step is alway 
trivial, it is most efficient to have the last sequence be as long as possible. This is why we order 
the sizes of the set partition by increasing order. 

Thus the sequence corresponding to the tree Tf^j is 

(7,8,2,3,2,4,3,3,2,1). 
It is easy to reconstruct ttx^ from this sequence and hence 

vfTi = ({8}, {10, 3}, {5, 4}, {9, 7}, {11, 6, 2}). 

It then follows that 

TTTi = (0, 0, {10, 3}, 0, {11, 6, 2}, 0, {5, 4}, 0, {8}, 0, {9, 7}, 0). 

The function /i corresponding to ttti and its image under Q are pictured in Figure 6. 

The problem of ranking and unranking trees with a given multiset of degrees is just an 
extension of the the problem ranking and unranking trees with a given sequence of degrees. 
That is, the distribution of degrees is just another set partition. For example, consider the 
sequence of degrees for the tree Tq pictured in Figure 5, s = (0, 0, 2, 0, 3, 0, 2, 0, 1, 0, 2, 0, 0). We 
can view that sequence as a set partition, A(s) = (Aq, . . . , A3) where Aj is the set of places 
where i appears in the sequence s. In our example, we would identify s with the set partition 

A(s) = ({1, 2, 4, 6, 8, 10, 12}, {9}, {3, 7, 11}, {5}). 

Just as in the case where we ranked and unranked set partitions associated with trees in 
Cn^g, it is more efficient if we rearrange the set partition by increasing size. This means that 
we must have a data structure to record the degrees associated with the set partition which in 
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Figure 6: The tree of rank 50,005 in Ci2,(o,o,2,o,3,o,2,o,i,o,2,o> 

this case is just the triples (Aj, |Ai|, i). It is easy to see that we can produce such a hst in Hnear 
time from the tree. In our example, we would produce the list 

A{s) = (({1, 2, 4, 6, 8, 10, 12}, 7, 0), ({9}, 1, 1), ({3, 7, 11}, 3, 2), ({5}, 1, 3)). 

Using a lexicographic bucket sort algorithm [Tl], we can sort this list according to the lexico- 
graphic order on the last two entries of the triples to produce the list 

AM = (({9}, 1, 1), ({5}, 1, 3), ({3, 7, 11}, 3, 2), ({1, 2, 4, 6, 8, 10, 12}, 7, 0)). 

Then we use this ordering to produce an ordered set partition vr^ where we ignore any empty 
partitions. In our example, we would produce 

= {9}, {5}, {3, 7, 11}, {1, 2, 4, 6, 8, 10, 12}). 

Finally, we use this set partition to produce a sequence of decreasing functions. That is, we let 
E = {1, . . . , n}. Let vf = {Ai, . . . ,Ak) be an ordered set partition of E where each block Ei 
is nonempty and ordered in decreasing order, Ai = {at^i > . . . > ai^a}, that comes from some 
element vr in vr^ as described above. Let (ri^i, . . . , ri^t^) be the ordered sequence of ranks of 
the respective (ai^i, . . . , ai^n) in E. In general, let (rj^i, . . . , rj^tj be the ranks of the respective 
(aj^i, . . . , tti^ti) in E \ U*~^ylj. For our example, 
(ri^i) = (9) can be considered an element of I'J^i2,i 
(^'2,1) = (5) can be considered an element of VFn^i 
(^3,1; '^3,2) ^^3,3) = (9, 6, 3) can be considered an element of I'-Fio.s 

(r4^i, r4_2) ^^4,3) ^^4,4) ^^4,5, '''4,6) ^"4,7) = (7, 6, 5, 4, 3, 2, 1) can be considered an element of VF'j^t^2- 
Thus we produce the sequence ifg = (9, 5, 9, 6, 3, , 7, 6, 5, 4, 3, 2, 1) to code the sequence s which 
can be considered an element of PJ^i2,i x x I'J^io,3 x VT-jj. We can then concatenate 

this sequence ifg with the sequence Tf^g to produce a sequence ^s.Tq ■ In our example, 

IlgTo = (9, 5, 9, 6, 3, , 7, 6, 5, 4, 3, 2, 16, 8, 5, 6, 2, 3, 1, 3, 2, 1) 

coming from an element of 

PJ"l2,l X VFii^i X PJ"l0,3 X VFtj X PJ"l0,l X PJ"9,2 X PJ"7,2 X PJ"5,2 X PJ"3,3 
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or as a leaf in the tree 

TvTu,l ® ^©^11,1 ® ^©^10,3 ® ^©^7,7 ® T-DT^o,! ® ^©^9,2 ® ^©^7,2 ® ^©^5,2 ® ^X'^3,3- 

Note that the size of the trees needed for the product lemma, Lemma [6j are 

l^2?^3,3l = (3) = ^' 

\Tvr,,,®TvT,,,\ = (2) • (3) =10' 

|ri,^7.2 ® TvT.a ® TvT,,,\ = (2) • (2) • (3) = 210, 

\TVT<,,2 ® ^25-^7,2 ® TVT5,2 ® ^15^3,3 ^©-^7,2 ® ^25-^5,2 ® ^©-^3,3 I = 

(^) • (D • • © = 7560, 

\TvTw 1 ^I'^g 2 ® T-Djr^ 2 Txjj-g 2 (g) rx,j-3 3 I = 
O-iD-Q-iD- it) =75600, ' 

|7i);f7 7 (g) Tx^j-jo 1 (g) Tpj-g 2 ® Tipj-^ 2 ^ ^©^5 2 ® ^25^3 3 1 = 
©•(T)©-©-©-© =75600, 

I^D^Fk, 3 (g) Tuj-^ 7 g) Tpj-^o ^ g) T^J-g 2 ® T-Djr^ 2 <^ ^©^5 2 ® ^I'^S 3 I = 

C3°)-Q-(\°)©-Q-(l(D =9,072,000; 

l^oj^ii 1 ^©j-^o 3 (g Tpj-^ 7 g) T-vjr^^ ^ g) Tx^j-g 2 <g Tdj-^ 2 (g Tbj-g 2 <g 7x,j-3 3 1 = 
(\')-C3°)-0-(\°)© = 99, 792, 000, 

|7i5;Fi2 1 ® ^©^ll 1 ® ^I'^lO 3 ® ^I'^7 7 <^ ^15^10 1 ^ '^■DT9 2 ® T-Djr^ 2 ^2?J-5 2 ® ^»^3 3 I = 

■ (V) • (^3°) • Q ■ (\°) • (D • il) ■ il) = 1, 197, 504, 000. 
Thus there are a total of 1,197,504,000 trees in Ci2,i whose degree sequence yields the multiset 
S = (0'', 1^, 2^, 3^). We can then use the product lemma. Lemma [6| to compute the rank of Tq 
in Ci2,s as follows. 

rank^^^ _^{To) = rankT)jr^^ -^{{9)) x 99,720,000 

+rankvru.A{^)) ^ 9,072,000 
+ranA;x)j-io.3((9,6,3)) x 75,600 
+ranfccj-7,7((7,6,5,4,3,2, 1)) x 75,600 
+ranA;x)j="io 1 ((5)) ^ 7,560 
+ran/cDj-9 2((8,5)) x 210 
+ranA;x)j7^ 2((6, 2)) x 10 
+ranA;Djrg 2((3, 1)) x 1. 

By Theorem [IJ we have that 
ranfci,^,2,i((9))= Ct') =8, 
rankvT,,,A{^)) = {\^) = 

rankvT^o,,{{9, 6, 3)) = {^f) + {\^) + (=^7^) = 56 + 10 + 2 = 68, 
ranfci,^, , ((7, 6, 5, 4, 3, 2,1)) = 0, 
rankvr,o,{{5)) = {''-^)=A, 

rankvT,,2 ((», 5)) = f + f 7^) == 21 + 4 = 25, 
rankvTr.2 ((6, 2)) = 1%^) + (^7^) = 10 + 1 = 11, 
rankT,^,,i{3, 1)) = (=^7^) + {\^) =1 + = 1. 
Thus 

rank{Wo) = (9 x 99, 720, 000) + (5 x 9, 072, 00) + (68 x 75, 600) + (0 x 75, 600) 
+(4 X 7560) + (25 x 210) + (11 x 10) + (1 x 1) = 843, 342, 641. 
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Thus the tree Tq pictured in Figure 5 has rank 843,342,641 among all the trees T G ^12, (0^,11, 23,31)- 
The unranking procedure for Cn,s comes from simply reversing the ranking procedure using 
Theorem [5] and Lemma [7| Again we will exhibit the procedure by finding the tree T2 whose 
rank is 60,000,00 in Ci2,(o'^,ii,2J',3i)- The first step is to carry out the series of quotients and 
remainder according to Lemma [7} In our case, this leads to the following calculations. 



60,000,00 = 


(6 X 99,792,000) + 1,248,000 


1,248,000 = 


(0 X 9,072,000) + 1,248,000 


1,248,000 = 


(16 X 75,600) + 38,400 


38,400 = 


(0 X 75,600) +38,400 


38, 400 = 


(5 X 7560) + 600 


600 = 


(2 X 210) + 180 


180 = 


(18 X 10) + 


= 


(0 X 1) + 0. 



It then follows that we can construct the sequence corresponding to vt^^^Tj by concatonating the 
sequences vi, . . . ,vq where 

1. vi is the decreasing function of rank 6 in 257^12, 1, 

2. V2 is the decreasing function of rank in DJ^n^i, 

3. V3 is the decreasing function of rank 16 in I'-Fics, 

4. V4 is the decreasing function of rank in VTjj and 

5. is the decreasing function of rank 5 in DJ-'io,i. 

6. ■yg is the decreasing function of rank 2 in T)Tg^2- 

7. v^ is the decreasing function of rank 18 in DT-j^i- 

8. Us is the decreasing function of rank in T>T^^2- 

9. "Ug is the decreasing function of rank in T>J^2,3- 

It is clear that the sequence of rank 6 in PJ-"i2,i is (7) and the sequence of rank in VJ^n^i is (1). 

To find the element (/(I), /(2), /(3)) of rank 16 in I'J^ics, we use the procedure in Theorem[5j 
We start by setting m' := 16 and s = 3. Since (3) = 10 < 16 < (3) = 20, then /(I) - 1 = 5 and 
hence /(I) = 6. Then we set m' := 16 — 10 = 1 and s = 2. Since (2) = 6 < 6 < (2) = 10, then 
/(2) — 1 = 4 and hence /(2) = 5. Finally we set m' := 6 — 6 = and ,5 = 1. Since — (5) = 0, 
we get that /(3) — 1 = or /(3) = 1. Thus (6, 5, 1) has rank 16 in DJ^io,3. 

It is clear that the sequence of rank in VTjj is (7, 6, 5, 4, 3, 2, 1). 

Next it is clear that element of rank 5 in VTiq^i is (6). 
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To find the element (/(l),/(2)) of rank 2 in PJ-"9,2, we use the procedure in Theorem [5| We 
start by setting m' := 2 and s = 2. Since (2) = 1 < ^ < (2) = 3, then /(I) — 1 = 2 and hence 
/(I) = 3. Then we set m' := 2 - 1 = 1 and s = 1. Since 1 - (J) = 0, we get that /(2) -1 = 1 
or /(2) = 2. Thus (3, 2) has rank 2 in 

To find the element (/(l),/(2)) of rank 18 in PJ^7_2, we again use the procedure in Theo- 
rem [sj We start by setting m! := 18 and s = 2. Since (2) = 15 < 18 < (2) = 21, then 
/(I) — 1 = 6 and hence /(I) = 7. Then we set m' := 18 — 15 = 3 and s = 1. Since 3 — Q = 0, 
we get that /(2) - 1 = 3 or /(2) = 4. Thus (7,4) has rank 18 in PJ"7,2- 

Finally the sequence of rank in DT-j^2 is clearly (2, 1) and the element of rank in DTz.i, 
is (3,2,1). 

Thus the sequences corresponding to the set partitions vf^ and tttt, are 

7f,- : (7, 1, 6, 5, 1, 7, 6, 5, 4, 3, 2, 1) and 
7fT2 : (6,3,2,7,4,2,1,3,2,1). 

It is easy to reconstruct vf^ and t^t^ to get that 

Tfs = ({7}, {1}, {8, 6, 2}, {12, 11, 10, 9, 5, 4, 3}) and 
TfTi = ({7}, {4, 3}, {11, 8}, {5, 2}, {10, 9, 6}) 

It then follows that 

s = (3,2,0,0,0,2,1,2,0,0,0,0) 

and 

TTT, = ({10, 9, 6}, {4, 3}, 0, 0, 0, {11, 8}, 0, {7}, 0, {5, 2}, 0, 0, 0). 

Thus the function corresponding to vrj-g and its image under are pictured in Figure 7. 

We note that essentially the same analysis of the complexity of ranking and unranking relative 
Cn,s applies to the complexity of ranking and unranking relative to Cn^s for any multiset S so 
that it requires 0(r?log(ri)) bit operations to rank and unrank relative Cn,s- 
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